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Preface 


The  purpose  of  this  report  is  to  describe  a  structured  method  that 
can  be  used  to  identify  the  potential  magnitude/s  and  significance 
of  meteorological  scale-decoupling  between  specific  surface  wind  ob¬ 
servation  sites  and  a  meso -ft  scale  prognostic  model.  This  work  is 
related  to  studies  that  have  examined  the  climatologically  optimal 
placement  of  surface  sensors  over  complex  terrain.  However,  results 
in  this  study  are  derived  using  a  meso-/3  scale  model  in  conjunction 
with  a  diverse  set  of  diurnally  varying,  nearly  meso-7  scale  surface 
wind  observation  data,  not  mean  climatological  data. 
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Executive  Summary 

Introduction 

Multiple  meso -/?  *  surface  wind  site  data  representativeness  is  exam¬ 
ined  here  in  a  comparison  of  corresponding  sets  of  surface  mesoscale 
model  windfield  predictions,  with  (experimental  case)  and  without 
(control  case)  periodic  surface  data  assimilation. 

Purpose 

This  investigation  was  conducted  to  identify  the  potential  magni¬ 
tude/s  and  significance  of  scale-decoupling  between  specific  observa¬ 
tion  sites  and  meso-/?  scale  model  prognostic  assumptions. 

Overview 

A  version  of  the  Higher  Order  Turbulence  Model  for  Atmospheric 
Circulation  (HOTMAC),  which  has  been  applied  to  a  large  variety  of 
meso-/?  scale  research  and  applications,  is  used  in  this  study  to  obtain 
selected  surface  gridded  windfields.  Six  different  radiosonde  initial¬ 
ization  data  sets  are  used  to  generate  HOTMAC  seven  hour  surface 
windfield  predictions  using:  1)  a  control  case  with  sparse  radiosonde 
initialization  data  (<  six  observations);  and  2)  an  experimental  case, 
in  which  supplementary  hourly  surface  wind  observations,  obtained 
from  a  consistent  set  of  15  non-uniformly  distributed,  automated  sur¬ 
face  observation  stations  over  a  meso-/?  scale  complex  terrain  area,  are 
also  periodically  assimilated  into  HOTMAC  windfields  by  locally  ad¬ 
justing  predicted  windfield  solutions  toward  observed  conditions  and 
correcting  for  mass  conservation  imbalances  in  the  adjusted  fields. 
Factor  separation  and  statistical  tests  are  used  to  characterize  the 
degree  of  surface  observation  station  scale-coupling  to  the  modeled 
meso-/?  scale  wind  field  forecasts. 

Conclusions 

The  mean  index  developed  for  use  in  this  study,  weighted  against 
testing  significance  levels ,  is  an  effective  means  with  which  extremely 
unrepresentative,  unpredictable  site  data  can  be  objectively  ranked, 
and  then  selectively  filtered  out  to  improve  the  dynamical  consistancy 
of  initialization  data  sets. 

*  Spatial  scale  ranges  are  approximately  2  km  <  Meso-{3  scale  <  200  km. 
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Recommendations 

The  diverse  range  of  individual  surface  observation  site  ensemble  com¬ 
parison  results  found  in  this  study  demonstrate  that  surface  data  sets 
require  careful  analysis  to  identify  each  site’s  scale-suitability  prior  to 
inclusion  in  a  meso -(3  scale  model  initialization  scheme. 


1.  Introduction 


Determining  optimal,  dynamically  consistent  initial  conditions  for 
making  predictions  is  the  goal  of  a  good  Numerical  Weather  Predic¬ 
tion  (NWP)  initialization  scheme.  However,  many  errors  can  be  intro¬ 
duced  into  initial  condition  data.  For  example,  data  void  regions  yield 
insufficient  data  density  to  resolve  mesoscale  features,  and  the  behav¬ 
ior  of  features  that  are  not  measured  will  probably  not  be  predicted 
very  well.  Other  possible  sources  of  initial  condition  errors  include 
observation  errors  resulting  from  instrument  performance  problems 
as  well  as  coding  and  transmission  mistakes.  The  representativeness 
of  observation  data  is  another  possible  source  of  initialization  condi¬ 
tion  errors. [1]  It  is  the  representativeness  problem  of  meso -0  scale  [2] 
initial  condition  and  periodic  assimilation  surface  wind  data  that  is 
investigated  in  this  study. 

This  study  centers  on  the  problem  of  meso -/?  surface  wind  data  rep¬ 
resentativeness.  Comparison  and  analysis  is  made  of  simulated  meso- 
j3  scale  anabatic  and  katabataic,  thermally  direct,  diurnally  induced 
and  gravity  forced  circulations  over  complex  terrain  with  and  without 
supplemental  surface  data  assimilation.  [3]  Meso-/3  scale  anabatic  and 
katabataic  circulations  that  interact  with  a  slowly  varying  ensemble 
flow,  with  otherwise  weak  synoptic  forcing  (i.e.,  the  main  scale  inter¬ 
action  t  in  this  study  is  essentially  limited  to  the  influence  of  mean 
flow  on  localized,  meso-/5  circulations). 

The  importance  of  mesoscale  interactions  was  recogonized  by  the  U.S. 

Army  as  early  as  1951  from  the  radar  research  results  of  Ligda  [2];  the 
term  mesoscale  was  actually  first  introduced  by  Ligda  in  a  U.S.  Army 
sponsored  work. [4]  The  U.S.  Army  today  has  a  continued,  increasing 
interest  in  understanding  more  about  mesoscale  interactions,  that  can 
have  a  significant  influence  on  battlefield  operations.  For  example, 
wind  variability  within  the  sphere  of  influence  of  the  mesoscale  can 
significantly  influence  artillery  accuracy. 

Accurate  short-range  windfield  forecasts  are  also  becoming  crucial  to 
other  U.S.  Army  field  operations,  such  as  in  predicting  atmospheric 
effects  on  chemical  fallout. [6]  As  a  result,  a  need  exists  for  the  U.S. 

Army  to  decrease  meteorological  uncertainty  and  errors  in  wind  field 

t  Identification  of  scale  differentiation  and  interactions  originated  primarily  from  comparing 
observations  of  atmospheric  phenomena  to  varying  observation  network  resolutions ;  this  was 
followed  by  making  theoretical  inferences  [4],  a  process  that  could  be  carried  out  using  factor 
separation. [5] 
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predictions.  Solutions  to  this  problem  translate  into  various  unit  per¬ 
formance  improvements:  for  example,  the  effectiveness  of  artillery 
fire.  [7]  Characterizing  the  representativeness  of  mesoscale  observa¬ 
tion  data  is  an  important  issue  (e.g.,  in  the  development  of  mesoscale 
windfield  prediction  decision  aids  for  the  U.S.  Army).  This  study  ex¬ 
amines  the  representativeness  problem  using  a  limited  set  of  surface 
wind  observation  data  that  is  selectively  incorporated  as  initial  con¬ 
dition  and  periodic  assimilation  data  into  a  mesoscale  atmospheric 
circulation  model. 

2.  Numerical  modeling  methods 


2.1  Description  of  the  model 

A  version  of  the  Higher  Order  Turbulence  Model  for  Atmospheric  Cir¬ 
culation  (HOTMAC),  a  second- moment  turbulence-closure  model,  is 
used  in  this  study  to  obtain  selected  4d  surface  windfield  gridded 
outputs.  The  HOTMAC  model  is  very  effective  in  predicting  meso- 
j3  scale  drainage  and  upslope  surface  winds  [8]  [9],  especially  under 
relatively  weak  synoptic  forcing;  under  these  conditions,  external  dif¬ 
ferential  surface  heating  effects,  which  are  estimated  by  HOTMAC, 
and  gravitational  forcing  over  complex  terrain  play  critical  roles. [3] 

HOTMAC  computes  geostrophic  wind  components  using  the  hydro¬ 
static  assumption  and  applies  a  Boussinesq  approximation  incorpo¬ 
rating  a  terrain  following  vertical  coordinate  system: 


z*  =  H 


^  Z  -  Zn 


H  —  Zg' 


(1) 


where  z *  and  z  are  the  transformed  and  Cartesian  vertical  coordi¬ 
nates,  respectively;  zg  represents  terrain  elevation  above  sea  level 
(ASL);  TI  is  the  material  surface  upper  boundary  of  the  model  in 
2*  coordinates  and  H  is  the  corresponding  z  coordinate  height. 


The  governing  conservation-of-momentum  equation  in  the  x  direction, 
following  the  coordinate  transformation,  is 


DU 

Dt 


f(V-Vg)  +  g!Lf- 

d  (  dU\  d_ 
dx  \  x  dx)  +  dy 


<  Oy  >\  8Zg 

6V  )  dx 


K xy 


H 

+  H-Zg 


(— uw ) 


—rj  Cd  [  a(z)  U  •  S  ]  +  G(Ut  —  U), 


(2) 
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where  <  >  indicates  an  average  over  a  horizontal  surface,  /  is  the 
Coriolis  parameter,  g  is  acceleration  due  to  gravity,  6V  represents 
virtual  potential  temperature  fluctuations,  and  Kx  and  Kxy  are  hori¬ 
zontal  eddy  viscosity  coefficients;  Cd  is  a  canopy  drag  coefficient,  a(z) 
is  the  plant  area  density,  rj  is  the  fraction  (0  <  r)  <  1)  of  area  covered 
by  trees,  and  S  —  \JU2  +  V2  +  W2.  The  pressure  gradient  force  term 
is  expressed  in  terms  of:  1)  Vg,  or  the  geostrophic  V  component;  and 
2)  the  second  term  on  the  right-hand  side  of  equation  (  2).  Mean 
turbulent  momentum  flux  is  represented  as  uw,  the  next  to  the  last 
term  on  the  right-hand  side  of  equation  (  2)  represents  drag  force  due 
to  tree  canopies,  and  the  last  term  nudges  the  horizontal  modeled 
wind  fields,  U,  to  target  winds,  [10]  * 

2.2  Case-specific  model  characteristics 

A  version  of  HOTMAC  was  configured  to  make  hourly  forecasts  out 
to  7  h  from  two  initialization  times  (0600L  and  1400L);  there  were 
18  x  22  modeled  horizontal  grid  points  and  31  model  layers  AGL  with 
an  upper  boundary  «  15  km  ASL.  The  southwest  terrain  anchor  was 
located  at  lat.  32.984  N,  long.  106.241  W,  with  8.0  km  model  grid 
resolution  both  in  the  x  and  y  directions.  Minimum  terrain  elevation 
was  1 ,185.0  m  ASL,  and  maximum  terrain  elevation  was  3,075  m  ASL. 

Only  available  radiosonde  data  were  used  to  initialize  HOTMAC  (see 
table  1).  Only  the  z  «  10  m  AGL  surface  windfield  outputs  were 
used  in  this  study  and  the  nudging  factor  G  was  held  constant  at 
G  =  5.0  x  10~4. 

2.3  Data  assimilation  modification  to  the  model 

A  modified  objective  analysis  procedure  was  used  to  convolve  grid- 
ded  HOTMAC  predictions  of  surface  windfields  with  corresponding, 
consistent  sets  of  randomly  spaced  surface  wind  observations.  This 
procedure  was  used  to  periodically  constrain  HOTMAC  windfield  pre¬ 
dictions  to  surface  wind  observation  data  as  it  became  available.  The 
goal  was  to  adjust  HOTMAC  windfield  predictions  (and  analyses  at 

t  Target  winds,  as  defined  in  [10],  are  obtained  by  solving  the  equations  of  motion;  for  ex¬ 
ample,  in  the  x  direction,  incorporating  observed  and  geostrophic  components  as  in 

Ut  =  Uob,  -  Jr(Vot,,  -  Vt),  (3) 

t-'n 

without  frictional  effects,  where  Cn  is  a  constant  nudging  coefficient.  Thus,  the  nudging  term 
adjusts  winds  toward  observation  data  in  the  free  atmosphere  but  has  little  effect  in  the  boundary 
layer  where  frictional  effects  become  significant.  Also,  Cn  =  G,  and  both  Cn  and  G  operate  in 
equation  (  2). 
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t0)  as  closely  as  possible  to  available  surface  wind  observations  with¬ 
out  eliminating  the  windfield  dynamics  predicted  by  HOTMAC  in 
data  sparse  regions.  This  procedure  was  used  to  perform  periodic 
reanalyses  of  hourly  HOTMAC  windfield  outputs. 


The  modified  objective  analyses  took  an  initial  HOTMAC  gridded 
surface  windfield  £  and  a  randomly  spaced  set  of  n  surface  wind  ob¬ 
servations.  [11]  The  mean  field  error  was  then  calculated  and  removed 
from  all  £  grid  points  to  produce  a  field  £*.  At  this  point,  all  £*  grid 
points  within  <  one  grid  square  were  adjusted  to  create  an  adjusted 
£**  field  defined  as  follows: 


r  =  r  + 


£"=1  Wm. 


(4) 


where  d*  is  the  local  discrepancy  between  observation  i  and  the  in¬ 
terpolated  £*  field  at  the  four  surrounding  local  grid  points  (m  = 
1,2, 3, 4);  Wmi  is  the  local  m’th  objective  weighting  factor.  Only  grid 
points  in  the  vicinity  of  the  randomly  spaced  surface  wind  observa¬ 
tions  are  affected  by  the  £**  operation.  The  unaffected  £*  grid  points 
are  then  adjusted  using  iterative  relaxation  to  alleviate  edge  disconti¬ 
nuity  between  C  and  £**  gridded  fields;  this  results  in  a  £***  gridded 
field  that  has  the  following  property 


v2r*  -  v2r, 


(5) 


The  purpose  of  this  objective  analyses  was  to  use  the  largest  set 
of  consistently  reported  surface  wind  observations  in  selectively  ad¬ 
justing  HOTMAC  windfield  predictions  to  corresponding  observation 
data. 


Reanalyzed  HOTMAC  surface  windfields  were  then  further  adjusted 
for  mass  consistency  by  iteratively  reducing  divergence,  where  diver¬ 
gence  is  defined  as  follows: 


du  dv  ...  .. 

o — b 
ox  oy 


(6) 


where  u  is  the  west-east  wind  component,  v  is  the  south-north  wind 
component,  and  <l>(i,j)  is  the  divergence  at  each  grid  point. [12]  The 
combination  of  the  modified  objective  analysis  and  mass  consistency 
operations  represent  the  data  assimilation  scheme  employed  in  this 
study.  Advantages  of  this  basic  assimilation  scheme  are  primarily  its 
economy  of  computational  resources,  processing  time,  and  input  data 
requirements. 
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3.  Meteorological  Scenarios 

Meteorological  observation  data,  consisting  of  radiosonde  and  auto¬ 
mated  surface  wind  observations  (~  10  m  AGL),  were  obtained  for 
three  separate  days  (Julian  days  167,  246,  and  272  in  1992)  in  a 
136  x  168  km2  area  centered  over  a  complex  terrain  region  (eleva¬ 
tion  extremes:  1,185.0  m  -  3,075.0  m  ASL;  absolute  terrain  differen¬ 
tial  =  1,890  m)  in  southern  New  Mexico  (see  figure  1).  This  data 
was  acquired  during  portions  of  a  Target  Area  Meteorological  Data 
Experiment  (TAMDE)  conducted  by  the  U.S.  Army  Research  Labo¬ 
ratory  (ARL),  Battlefield  Environment  Directorate  (BED).  TAMDE 
was  implemented  to  study  the  variability  of  weather  and  weather  ef¬ 
fects  on  a  variety  of  weapon  and  fire  control  systems.  [6] 

TAMDE  radiosonde  data  used  to  initialize  the  control  and  experi¬ 
mental  HOTMAC  simulations  are  identified  in  table  1.  Three  sets  of 
hourly  surface  wind  observation  data,  corresponding  to  Julian  days 
167,  246,  and  272  in  1992,  were  collected  for  the  experimental  assimi¬ 
lation  portion  of  this  study.  The  10  m  AGL  surface  wind  observations 
were  obtained  from  a  consistent  network  of  15  automated  surface  ob¬ 
servation  sites  (relative  sensor  locations  indicated  in  figure  1)  that 
were  configured  to  log  15  minute  averaged  observation  data. 

4.  Experiment  Design 

4.1  Control  case 


In  the  control  case,  hourly  surface  wind  data  were  objectively  ana¬ 
lyzed  using  the  following  basic  interpolation  scheme: 

U{i,j)  =  ^ m(r)Um/  Mr),  (7) 

771  =  1  771=1 


where  U(i,j )  is  the  calculated  value  at  the  grid  point  (i,j),  Um  is 
the  observed  wind  at  site  m,  ^m(r)  =  ^ ,  and  r  is  the  grid  point  to 
observation  site  distance .  Following  the  1  /r2  interpolation  of  observed 
surface  winds  to  gridded  fields,  the  gridded  fields  were  adjusted  for 
mass  consistency  using  iteratively  reduced  divergence  to  a  predefined 
threshold,  where  divergence  is  defined  as  follows: 


du  dv  ..  .. 

ai  +  Ty=^3)' 


(8) 


where  u  is  the  west-east  wind  component,  v  is  the  south-north  wind 
component,  and  (f>(i,j)  is  the  divergence  at  each  grid  point. [12] 
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4.2  Experimental  case 

In  the  experimental  case  a  modified  objective  analysis  procedure,  sum¬ 
marized  in  equations  (  4-  6),  was  used  to  convolve  gridded  HOTMAC 
surface  windfield  predictions  with  corresponding  sets  of  irregularly 
spaced  surface  wind  observations;  this  periodically  constrained  HOT¬ 
MAC  windfield  predictions  to  surface  wind  observation  data.  There¬ 
fore,  resulting  surface  windfield  predictions  were  a  function  of  both:  1 ) 
HOTMAC  meso -/?  prognostic  atmospheric  dynamics  considerations; 
and  2)  the  selected  surface  observation  data. 

4.3  Data  analysis  techniques 

Because  anabatic  (katabatic)  circulation  interactions  over  a  complex 
terrain  domain  often  produce  confluent  (diffluent)  flow  in  valleys 
(along  ridgelines),  comparing  individual  points  in  HOTMAC  surface 
windfield  predictions  to  mean  windfield  statistical  properties  is  of 
little  assistance  in  characterizing  the  true  meso-/?  scale  representa¬ 
tiveness  of  a  particular  observation  site.  Instead,  the  corresponding 
control  and  experimental  case  results  need  to  be  compared  at  each 
point  that  contributed  observation  data  to  this  study. 

Factor  separation  was  used  to  compute  interactions  of  the  experimen¬ 
tal  case  factor  influencing  HOTMAC  surface  windfield  predictions.  [5] 
The  interaction  factor  investigated  in  this  study  is  meso-/?  scale  sur¬ 
face  wind  unpredictability,  A,  which  is  inversely  related  to  observation 
site  representativeness;  A  is  the  effect  of  periodically  assimilated  sur¬ 
face  data  on  a  meso-/?  scale  HOTMAC  base  case  simulation  field,  /( 0). 
Factor  A  is  multiplied  by  a  changing  coefficient  c,  where  0  <  c  <  1  in 
the  following: 

A(c)  =  cA,  (9) 

which  produces  a  resulting  field  /  that  is  a  continuous  function  of  c: 

/  =  /(c),  (10) 

where  function  /(c)  can  be  decomposed  into:  1)  a  constant  part  that 
is  independent  of  c,  /0,  which  is  a  short  form  of  /( 0)  and  represents  the 
fraction  of  /  that  does  not  depend  on  factor  A;  and  2)  a  c  dependent 
component,  /(c). 

In  this  case,  the  surface  windfield  unpredictability,  /(c),  can  be  ex¬ 
pressed  as  follows: 

m  =  /(c)  -  (ii) 

where  /(c)  is  the  value  of  a  predicted  experimental  case  HOTMAC 
surface  windfield  that  is  a  continuous  function  of  c,  and  f0  is  the  value 
of  /  in  the  control  simulation  that  has  factor  A  omitted. 
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Table  1.  TAMDE  radiosonde  data  used  to  initialize  HOTMAC  con¬ 
trol  and  experimental  simulations,  (availability  indicated  with  a  • 
symbol). 


UA  Station  No. 

DATE 

TIME 

1 

2 

3 

4 

5 

6 

167 

0600 

• 

1400 

• 

246 

0600 

• 

• 

• 

1400 

• 

• 

• 

• 

• 

272 

0600 

• 

• 

• 

• 

1400 

• 

• 

• 

• 

• 

• 

Figure  1 .  Model  surface  terrain  elevation  ASL  (contour  interval  =  250  m, 
area  =  136  km  x  168  km).  Automated  surface  (radiosonde)  obser¬ 
vation  site  locations  are  indicated  by  corresponding  site  alpha  (*  nu¬ 
meric)  symbols. 
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In  this  study,  /(c)  is  statistically  analyzed  to  isolate  and  resolve  the 
inversely  related  representativeness  of  selected  observation  site  data 
to  the  simulated  HOTMAC  meso -/?  scale  regime  of  anabatic/katabatic 
boundary  layer  circulation  interactions  using  the  following  measures 
of  the  /(c)  surface  u  and  v  wind  component:  1)  correlation  and  scat¬ 
ter;  2)  relative  mean  bias;  3)  percentage  of  cases  within  a  factor  of 
two;  4)  mean  and  upper/lower  bound  differences;  5)  F-test  for  sig¬ 
nificance  of  variance  differences;  6)  Student’s  t-test  for  significance  of 
differences  in  means;  7)  Kolmogorov-Smirnov  test  for  significance  in 
cumulative  distribution  function  differences;  and  8)  Chi-Square  test 
for  significance  in  population  distribution  differences.  [13]  [14]  [15] 


5.  Results 

The  /(c)  =  \J{f{c)u\2  +  [ f(c)v ]2  combined  results  for  each  surface 
observation  site  (see  figure  1)  are  illustrated  in  figure  2,  where  R  = 
/(c),  and  the  following  is  true: 


f(c)u  —  ^ experimental  ■  ^control 5 
/(c)v  —  V  experimental  ^ control • 


(12) 

(13) 


Results  from  this  test  indicate  that  sites  1,4,  and  11  had  the  smallest 
/(c)  differences  (<  2.4  m/s)  and  are,  therefore,  the  most  representa¬ 
tive  meso -/?  scale  sites.  Conversely,  sites  8,  10,  and  13  had  the  largest 
/(c)  differences  (>  3.0  m/s),  and,  according  to  this  test  result,  were 
the  least  representative  meso-/?  scale  sites. 


The  R  comparisons  in  figure  2  were  integrated  into  the  measure  of 
/(c)  mean  and  upper/lower  bound  differences  (test  number  4,  from 
section  4.c).  The  remaining  seven  f(c)u  and  /(c)„  tests  mentioned 
in  section  4.c  were  applied  to  /(c)  data  interpolated  at  each  surface 
observation  site  (figure  1)  from  control  and  experimental  simulation 
results.  All  significant  test  results  are  summarized  in  table  2,  where 
a  (“(g)”)  symbol  indicates  high  confidence  in  identifying  relatively 
small  (large)  values  of  /(c);  small,  relative  /(c)  results  for  particular 
tests,  at  particular  sites,  indicate  a  higher  degree  of  meso-/?  scale 
representativeness. 

A  simple  mean  index  scheme  was  devised  to  summarize  the  net  test 
results  for  each  observation  site  (see  table  2).  The  mean  index  for 
each  site  was  calculated  by  dividing  the  number  of  representative 
test  results  (“•”  symbols  from  table  2)  by  the  sum  of  all  high  confi¬ 
dence  test  result  relative  rankings  (“•”  and  “(g)”  symbols  from  table 
2).  Therefore,  higher  (lower)  mean  index  values  are  associated  with 
relatively  more  representative  (unrepresentative)  sites.  For  example, 
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sites  14  and  15  have  perfect  mean  index  values  (1.0)  and,  without 
considering  any  other  factors,  are  the  most  representative  observa¬ 
tion  sites.  Conversely,  sites  3  and  12  had  mean  in  dear  values  <  0.17, 
indicating  the  lowest  degree  of  relative  meso -/?  scale  observation  site 
agreement. 

However,  mean  index  values  are  also  a  function  of  the  total  number 
of  high  confidence  test  results  per  observation  site.  Thus,  the  mean 
index  values  require  further  interpretation.  For  example,  the  mean 
index  can  be  weighted  against  its  significance  level.  For  the  purposes 
of  this  study,  a  mean  index  -  significance  level  was  defined  as  the 
cumulative  number  of  high  confidence  test  results  divided  by  the  total 
number  of  tests  used  (n  =  8).  Table  2  lists  significance  levels  and 
mean  index  values  for  each  site,  and  figure  3  presents  a  graphical 
comparison  of  both  significance  level  and  the  mean  index  for  each 
site.  Figure  3  indicates  that  site  3  scored  the  highest  significance 
level  (1.0),  while  sites  5,  14,  and  15  had  the  lowest  significance  levels 
(0.25). 

6.  Summary  and  Conclusions 

This  study  was  implemented  to  investigate  the  representativeness  of 
operationally  available  data.  The  goal  was  to  identify  how  well  avail¬ 
able  observations  from  a  given  set  of  observation  sites  compared  to 
reliably  predicted  meso-/?  scale  circulation  features.  The  mean  index 
derived  in  this  study  (see  table  2)  is  a  simple  approach  to  the  com¬ 
plex  quality  control  problem  of  classifying  relative  representativeness 
of  meso-/?  scale  observation  sites  in  a  complex  terrain  setting.  The 
mean  index,  defined  in  section  5,  can  be  used  to  characterize  the 
spectrum  of  individual  observation  site  agreement  with  dynamicaly 
balanced  katabatic/anabatic  meso-/?  scale  modeled  circulation  inter¬ 
action  in  complex  terrain  settings;  this  provides  objective  guidance 
on  a  particular  observation  site’s  meso-/?  scale  representativeness. 

The  /(c)  results  in  figure  2  agreed  very  well  with  the  mean  index 
results  (see  table  2).  For  example,  sites  1,  4,  and  11  had  the  lowest  R 
values  (figure  2)  and  all  three  of  these  sites  scored  significantly  above 
the  average  mean  index  results  (.83,  .75,  and  .75,  respectively,  with 

>  4  independent  test  indicators,  or  having  significance  levels  >  .50); 
these  combined  results  characterize  sites  1,  4,  and  11  as  the  most 
representative  at  the  meso-/?  scale.  Conversely,  sites  3,  8,  and  13 
had  below  average  mean  index  results  (.13,  .33,  and  .43,  respectively, 
with  >  6  independent  test  indicators,  or  having  significance  levels 

>  .75),  and  all  three  of  these  sites  had  the  highest  R  values  (figure 

2). 
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For  site  10,  R  >  3.0  m/s,  which  was  greater  than  R  results  for  sites 
3  and  12  (see  figure  2).  However,  only  three  out  of  eight  possible 
independent  test  indicators  in  table  2  could  clearly  differentiate  site 
10  relative  representativeness  ( significance  level  of  .38).  Therefore, 
site  10  results  were  not  considered  in  further  relative  site  rankings. 
Site  12  did  not  have  a  greater  R  value  than  sites  3,  8,  or  13  (figure 
2);  in  contrast  to  site  10  inferences,  site  12  scored  the  second  lowest 
mean  index  (.17  —  determined  by  6  test  indicators,  a  signifcance  level 
of  .75),  and  was  classified  as  an  unrepresentative  observation  site. 
Thus,  the  least  representative  observation  sites  in  this  study  include 
sites  3,  8,  12,  and  13.  These  results  illustrate:  1)  why  site  representa¬ 
tiveness  conclusions  should  not  be  dependent  on  only  one  test  result; 
2)  how  the  cumulative  comparisons  of  the  number  of  test  indicators 
can  become  a  significant  analysis/ decision  factor. 

From  figure  1,  sites  12  (1)  and  13  (m)  are  located  along  the  east¬ 
ern  slopes  of  mountainous  ridgelines  with  terrain  elevations  varying 
>  500  m  within  a  1  km  radius  of  each  of  these  sites.  Unpredictable 
mechanical  turbulance  effects,  forced  by  the  significant  variations  in 
local  terrain  elevation,  are  probable  causes  for  the  relatively  unrepre- 
sententative  ranking  of  these  two  sites.  Similarly,  site  3  (c)  is  within 
«  2  km  of  steep  terrain  rise  (>  750  m)  to  the  west.  Site  8  (h)  is  a 
difficult  observation  site  in  this  study  because  of  its  vicinity  to  the 
northern  model  boundary  condition  (<  1  km  within  the  model  do¬ 
main);  in  this  case,  the  observation  site  may  be  more  accurate  at 
the  meso-fi  scale  than  simulation  results  approaching  model  domain 
limits. 

Conversely,  sites  4  (d)  and  11  (k)  are  nested  >  3  km  within  the 
model  domain,  and  are  located  on  relatively  smooth  local  terrain 
with  minimal  relative  slope  changes;  this  is  probably  why  these  two 
sites  ranked  as  highly  as  representative  at  the  meso-fi  scale  in  terms  of 
their  corresponding  R,  mean  index ,  and  significance  level  Similar  to 
site  3  (c),  site  1  (a)  is  within  ft:  2  km  of  a  steep  change  in  terrain  (>  1 
km),  yet  like  sites  4  (d)  and  11  (k),  it  is  located  on  relatively  smooth 
local  terrain,  and  also  ranked  as  highly  representative  at  the  meso-/3 
scale;  Site  1  (a)  thus  illustrates  why  spatial  terrain  analysis  alone 
provides  insufficient  guidance  on  the  meso-fi  scale  boundary  layer 
representativeness  of  particular  observation  sites  when  other  analysis 
tools,  like  a  valid  meso-fi  scale  atmospheric  circulation  model,  are 
available. 


Table  2.  Site/test  summary  matrix,  identifying  the  most  representa¬ 
tive  (least  representative)  meso -/?  scale  observation  sites  with  a 
(“®”)  symbol,  where  test  numbers  correspond  to  significant  u  and/or 
v  surface  wind  component  /(c)  measures.^  The  significance  level  and 
mean  index  are  defined  in  section  5. 


SITE 

2 

n 

1 

11 

6 

7 

8 

SIGNIFICANCE 

LEVEL 

MEAN 

INDEX 

1 

4-* 

a 

• 

• 

• 

• 

a 

.75 

.83 

2 

4-> 

b 

• 

.63 

3 

4-> 

c 

El 

• 

1 

1.0 

.13 

4 

4— > 

d 

• 

• 

• 

.50 

.75 

5 

4— ► 

e 

• 

.25 

.50 

6 

4— > 

f 

• 

EM 

• 

.50 

.50 

7 

4— ► 

g 

• 

• 

.38 

.67 

8 

4— y 

a 

• 

EM 

E 

• 

.75 

.33 

9 

4 — > 

i 

• 

• 

• 

.50 

.75 

10 

4-* 

j 

0 

EM 

> 

.38 

.00 

11 

4— > 

k 

• 

D 

• 

E 

.50 

.75 

12 

4-)- 

11 

151 

0 

• 

.75 

.17 

13 

4— > 

m 

D 

EM 

El 

n 

I 

D 

.88 

.43 

14 

4 — ► 

n 

• 

• 

.25 

1.00 

15 

4 — y 

o 

• 

• 

.25 

1.00 

£  USAGE 

9 

D 

7 

10 

6 

JLJ 

8 

AVG=0.54 

f  NOTE:  Test  Numbers  correspond  to  the  following  statistical  methods: 

1  -  correlation  and  scatter; 

2  -  relative  mean  bias; 

3  -  percentage  of  cases  within  a  factor  of  two; 

4  -  mean  and  upper/lower  bound  differences; 

5  -  F-test  for  significance  of  variance  differences; 

6  -  Student’s  t-test  for  significance  of  differences  in  means; 

7  -  Kolmogorov- Smirnov  test  for  significance  in  cumulative  distribution 

function  differences; 

8  -  Chi-Square  test  for  significance  in  population  distribution  differences. 
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SITE  NUMBER 


Figure  2.  Mean  R  values,  defined  in  section  5,  with  ±1.0  standard  deviation 
tick  marks,  for  each  surface  wind  observation  site.  Mean  R  for  all  sites  is  also 
included  (horizontal  dashed  line). 


Figure  3.  Mean  index  (*)  and  significance  level  (bar  chart  plot)  for  each  site, 
with  the  net  average  value  of  each  (plotted  respectively  as  a  dashed  and  dotted 
horizontal  line). 
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The  diverse  range  of  individual  surface  observation  site  ensemble  com¬ 
parison  results  in  this  study  demonstrate  that  surface  data  sets  re¬ 
quire  careful  analysis  to  identify  each  site’s  representativeness  prior 
to  inclusion  in  a  meso -/3  scale  model  initialization  scheme.  The  mean 
index  developed  for  use  in  this  study,  weighted  against  testing  signifi¬ 
cance  levels,  is  an  effective  means  with  which  extremely  unrepresenta¬ 
tive,  unpredictable  site  data  can  be  objectively  ranked;  physically  un¬ 
reasonable  site  data  can  then  be  selectively  filtered  out  to  improve  the 
dynamical  consistancy  of  initialization  data  sets,  which  will  decrease 
prediction  errors  when  a  meso  -ft  scale  model  is  integrated  forward  in 
time.[l][16] 
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Recommendations 


Critical  field  application/s  of  meso-(3  scale  model  forecasts  should  be 
accompanied  by  some  type  of  initialization  data  quality  control  that 
addresses  the  scale-coupling  of  observation  data  to  the  model  assump¬ 
tions;  otherwise,  forecast  errors  may  result  that  are  not  due  to  the 
model.  Further  classification  of  observation  site  data  should  also  in¬ 
clude  the  application  of  more  sophisticated  complex  quality  control 
schemes  that  more  specifically  address  other  factors  such  as  instru¬ 
ment  errror  and/or  calibration  problems;  these  will  be  implemented 
in  the  planned  continuation  of  this  research  work.§ 


§  Brock  [14]  suggests  some  methodology  to  identify  suspect  observation  data  in  a  preliminary 
quality  assurance  audit  prior  to  its  incorporation  into  a  model  initialization  scheme.  1)  detailed 
observation  variability  intercomparisons;  2)  data  range  testing  against  preset  parameter  limits; 
and  3)  applying  an  impulse  noise  threshold  detection  filter  to  the  data. 


References 

1.  Carr,  F.,  “Introduction  to  Numerical  Weather  Prediction  and  Present/Future 
Modeling  Activities  at  NMC,”  Lecture  Notes  ( November  2,  1993),  COMET 
Mesoscaie  Forecasting  and  Analysis  Course,  NCAR,  Boulder,  CO  80307,  70 
pp.,  1993. 

2.  Fujita,  T.T.,  “Mesoscaie  Classifications:  Their  History  and  Their  Appli¬ 
cation  to  Forecasting,”  Mesoscaie  Meteorology  and  Forecasting ,  Ray,  P.S., 

Ed.,  Amer.  Meteor.  Soc.,  18  pp.,  1986. 

3.  Huschke,  R.,  Ed.,  “Glossary  of  Meteorology.”  Amer.  Meteor.  Soc.,  638 
pp.,  1959. 

4.  Emanuel,  K.A.,  “Overview  and  Definition  of  Mesoscaie  Meteorology,”  Mesoscaie 
Meteorology  and  Forecasting ,  Ray,  P.S.,  Ed.,  Amer.  Meteor.  Soc.,  17  pp., 
1986. 

5.  Stien,  U.,  and  Alpert,  P.,  “Factor  Separation  in  Numerical  Simulations,” 

J.  Atmos.  Sci.,  50,  9  pp.,  1993. 

6.  Grace,  J.,  “TAMDE  -  The  Variability  of  Weather  over  an  Army  Division 
Size  Area,”  1993  Battlefield  Atmospherics  Conference  Proceedings  (under 
final  revision),  Battlefield  Environment  Directorate,  White  Sands  Missile 
Range,  NM  88002,  13  pp.,  1993. 

7.  Blanco,  A.,  and  Cogan,  J.,  “Mobile  Profiler  System  Improvements  to  the 
MET  Error  Contribution  of  the  Artillery  Error  Budget,”  Battlefield  At¬ 
mospherics  Conference  Proceedings,  Battlefield  Environment  Directorate, 
White  Sands  Missile  Range,  NM  88002,  15  pp.,  1993. 

8.  Yamada,  T.,  “A  Numerical  Model  Study  of  Turbulent  Airflow  in  and  Above 
a  Forest  Canopy,”  J.  Meteor.  Soc.  of  Japan,  60,  15  pp.,  1982. 

9.  Yamada,  T.,  and  Bunker,  S.,  “A  Numerical  Model  Study  of  Nocturnal 
Drainage  Flows  with  Strong  Wind  and  Temperature  Gradients,”  J.  Appl. 
Meteor.,  28,  9  pp.,  1989. 

10.  Henmi,  T.,  Dumais,  R.,  and  Smith,  T,  “Operational  Short-Range  Forecast 
Model  for  Battlescale  Areas,”  Battlefield  Atmospherics  Conference  Pro¬ 
ceedings,  Battlefield  Environment  Directorate,  White  Sands  Missile  Range, 
NM  88002,  10  pp.,  1993. 

11.  Danard,  M.,  Lyv,  G.,  and  MacGillivray,  G.,  “A  Mesoscaie  Bulk  Model 
of  the  Atmospheric  Boundary  Layer,  Part  I,  Phases  I  and  II,”  Technical 
Report  ASL-CR-85-0126-1,  Contract  No.  DAAD07-83-C-0126,  U.S.  Army 
Atmospheric  Sciences  Laboratory,  White  Sands  Missile  Range,  NM  88002, 

61  pp.,  1985. 


21 


12.  Lee,  M.,  “Complex  Terrain  Wind  Model  Evaluation,”  Thesis ,  New  Mexico 
State  University,  Las  Cruces,  NM  88003,  149  pp.,  1992. 

13.  Hanna,  S.,  “Mesoscale  Meteorological  Model  Evaluation  Techniques,  with 
Emphisis  on  Needs  of  Air  Quality  Models,”  contained  in  AMS  Monograph 
(under  final  revision )  on  the  1992  Mesoscale  Workshop,  El  Paso,  TX,  19 
pp.,  1992. 

14.  Brock,  F.  V.,  “Analysis  of  the  Surface  Atmospheric  Measurement  System, 
Scientific  Services  Program  Report,”  Contract  No.  DAAL03-86-D-001, 
U.S.  Army  Atmospheric  Sciences  Laboratory,  White  Sands  Missile  Range, 
NM  88002,  50  pp.,  1990. 

15.  Press,  W.,  Flannery,  B.,  Teukolsky,  S.,  and  Vetterling,  W.,  “Numerical 
Recipes  in  C,”  Cambridge  University  Press,  735  pp.,  1988. 

16.  Toth,  Z.,  and  Kalnay,  E.,  “Ensemble  Forecasting  at  NMC:  The  Generation 
of  Perturbations,”  Bull.  Amer.  Meteor.  Soc .,  74,  14  pp.,  1993. 


Acronyms  and  Abreviations 


AGL 

ARL 

ASL 

BED 

HOTMAC 

NWP 

TAMDE 


Above  Ground  Level 
Army  Research  Laboratory 
Above  Sea  Level 

Battlefield  Environment  Directorate 

Higher  Order  Turbulence  Model  for  Atmospheric  Circulation 

Numerical  Weather  Prediction 

Target  Area  Meteorological  Data  Experiment 
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Kenneth  R.  Hardy 

ORG/91-Ol  B/255  1 

3251  Hanover  Street 
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